Genome-wide association study identifies new HLA class II haplotypes strongly protective against narcolepsy

Narcolepsy is a rare sleep disorder with the strongest human leukocyte antigen (HLA) association ever reported. Since the associated HLA-DRB1*1501-DQB1*0602 haplotype is common in the general population (15-25%), it has been suggested that it is almost necessary but not sufficient for developing narcolepsy. To further define the genetic basis of narcolepsy risk, we performed a genome-wide association study (GWAS) in 562 European individuals with narcolepsy (cases) and 702 ethnically matched controls, with independent replication in 370 cases and 495 controls, all heterozygous for DRB1*1501-DQB1*0602. We found association with a protective variant near HLA-DQA2 (rs2858884; P < 3 x 10(-8)). Further analysis revealed that rs2858884 is strongly linked to DRB1*03-DQB1*02 (P < 4 x 10(-43)) and DRB1*1301-DQB1*0603 (P < 3 x 10(-7)). Cases almost never carried a trans DRB1*1301-DQB1*0603 haplotype (odds ratio = 0.02; P < 6 x 10(-14)). This unexpected protective HLA haplotype suggests a virtually causal involvement of the HLA region in narcolepsy susceptibility.

Common CFTR haplotypes and susceptibility to chronic pancreatitis and congenital bilateral absence of the vas deferens

CFTR mutations enhance susceptibility for idiopathic chronic pancreatitis (ICP) and congenital bilateral absence of the vas deferens (CBAVD); however, it is unknown why CFTR heterozygotes are at increased disease risk. We recently showed that common CFTR variants are associated with aberrantly spliced transcripts. Here, we genotyped for common CFTR variants and tested for associations in two ICP (ICP-A: 126 patients, 319 controls; ICP-B: 666 patients, 1,181 controls) and a CBAVD population (305 patients, 319 controls). Haplotype H10 (TG11-T7-470V) conferred protection (ICP-A: OR 0.19, P<0.0001; ICP-B: OR 0.78, P = 0.06; CBAVD OR 0.08, P<0.001), whereas haplotype H3 (TG10-T7-470M) increased disease risk (ICP-A: OR 8.34, P = 0.003; ICP-B: OR 1.88, P = 0.007; CBAVD: OR 5.67, P = 0.01). The risk of heterozygous CFTR mutations carriers for ICP (OR 2.44, P<0.001) and CBAVD (OR 14.73, P<0.001) was fully abrogated by the H10/H10 genotype. Similarly, ICP risk of heterozygous p.Asn34Ser SPINK1 mutation carriers (OR 10.34, P<0.001) was compensated by H10/H10. Thus, common CFTR haplotypes modulate ICP and CBAVD susceptibility alone and in heterozygous CFTR and p.Asn34Ser mutation carriers. Determination of these haplotypes helps to stratify carriers into high- and low-risk subjects, providing helpful information for genetic counseling.

Globally dispersed Y chromosomal haplotypes in wild and domestic sheep

To date, investigations of genetic diversity and the origins of domestication in sheep have utilised autosomal microsatellites and variation in the mitochondrial genome. We present the first analysis of both domestic and wild sheep using genetic markers residing on the ovine Y chromosome. Analysis of a single nucleotide polymorphism (oY1) in the SRY promoter region revealed that allele A-oY1 was present in all wild bighorn sheep (Ovis canadensis), two subspecies of thinhorn sheep (Ovis dalli), European Mouflon (Ovis musimon) and the Barbary (Ammontragis lervia). A-oY1 also had the highest frequency (71.4%) within 458 domestic sheep drawn from 65 breeds sampled from Africa, Asia, Australia, the Caribbean, Europe, the Middle East and Central Asia. Sequence analysis of a second locus, microsatellite SRYM18, revealed a compound repeat array displaying fixed differences, which identified bighorn and thinhorn sheep as distinct from the European Mouflon and domestic animals. Combined genotypic data identified 11 male-specific haplotypes that represented at least two separate lineages. Investigation of the geographical distribution of each haplotype revealed that one (H6) was both very common and widespread in the global sample of domestic breeds. The remaining haplotypes each displayed more restricted and informative distributions. For example, H5 was likely founded following the domestication of European breeds and was used to trace the recent transportation of animals to both the Caribbean and Australia. A high rate of Y chromosomal dispersal appears to have taken place during the development of domestic sheep as only 12.9% of the total observed variation was partitioned between major geographical regions.

Influence of ABCB1, ABCC1, ABCC2, and ABCG2 haplotypes on the cellular exposure of nelfinavir in vivo

OBJECTIVES: The human immunodeficiency virus protease inhibitor nelfinavir is substrate of polyspecific drug transporters encoded by ABCB1 (P-glycoprotein), ABCC1 (MRP1) and ABCC2 (MRP2), and an inhibitor of BCRP, encoded by ABCG2. Genetic polymorphism in these genes may be associated with changes in transport function. METHODS: A comprehensive evaluation of single nucleotide polymorphisms (39 SNPs in ABCB1, 7 in ABCC1, 27 in ABCC2, and 16 in ABCG2), and inferred haplotypes was done to assess possible associations of genetic variants with cellular exposure of nelfinavir in vivo. Analysis used peripheral mononuclear cells from individuals receiving nelfinavir (n=28). Key results were re-examined in a larger sample size (n=129) contributing data on plasma drug levels. RESULTS AND CONCLUSIONS: There was no significant association between cellular nelfinavir area under the curve (AUC) and SNPs or haplotypes at ABCC1, ABCC2, ABCG2. There was an association with cellular exposure for two loci in strong linkage disequilibrium: ABCB1 3435C>T; AUCTT>AUCCT>AUCCC (ratio 2.1, 1.4, 1, Ptrend=0.01), and intron 26 +80T>C; AUCCC> AUCCT > AUCTT (ratio 2.4, 1.3, 1, Ptrend=0.006). Haplotypic analysis using tagging SNPs did not improve the single SNP association values.

Accumulating mutations in series of haplotypes at the KIT and MITF loci are major determinants of white markings in Franches-Montagnes horses.

Coat color and pattern variations in domestic animals are frequently inherited as simple monogenic traits, but a number are known to have a complex genetic basis. While the analysis of complex trait data remains a challenge in all species, we can use the reduced haplotypic diversity in domestic animal populations to gain insight into the genomic interactions underlying complex phenotypes. White face and leg markings are examples of complex traits in horses where little is known of the underlying genetics. In this study, Franches-Montagnes (FM) horses were scored for the occurrence of white facial and leg markings using a standardized scoring system. A genome-wide association study (GWAS) was performed for several white patterning traits in 1,077 FM horses. Seven quantitative trait loci (QTL) affecting the white marking score with p-values p≤10(-4) were identified. Three loci, MC1R and the known white spotting genes, KIT and MITF, were identified as the major loci underlying the extent of white patterning in this breed. Together, the seven loci explain 54% of the genetic variance in total white marking score, while MITF and KIT alone account for 26%. Although MITF and KIT are the major loci controlling white patterning, their influence varies according to the basic coat color of the horse and the specific body location of the white patterning. Fine mapping across the MITF and KIT loci was used to characterize haplotypes present. Phylogenetic relationships among haplotypes were calculated to assess their selective and evolutionary influences on the extent of white patterning. This novel approach shows that KIT and MITF act in an additive manner and that accumulating mutations at these loci progressively increase the extent of white markings.

Multiple regulatory variants located in cell type-specific enhancers within the PKP2 locus form major risk and protective haplotypes for canine atopic dermatitis in German shepherd dogs.

BACKGROUND Canine atopic dermatitis (CAD) is a chronic inflammatory skin disease triggered by allergic reactions involving IgE antibodies directed towards environmental allergens. We previously identified a ~1.5 Mb locus on canine chromosome 27 associated with CAD in German shepherd dogs (GSDs). Fine-mapping indicated association closest to the PKP2 gene encoding plakophilin 2. RESULTS Additional genotyping and association analyses in GSDs combined with control dogs from five breeds with low-risk for CAD revealed the top SNP 27:19,086,778 (p = 1.4 × 10(-7)) and a rare ~48 kb risk haplotype overlapping the PKP2 gene and shared only with other high-risk CAD breeds. We selected altogether nine SNPs (four top-associated in GSDs and five within the ~48 kb risk haplotype) that spanned ~280 kb forming one risk haplotype carried by 35 % of the GSD cases and 10 % of the GSD controls (OR = 5.1, p = 5.9 × 10(-5)), and another haplotype present in 85 % of the GSD cases and 98 % of the GSD controls and conferring a protective effect against CAD in GSDs (OR = 0.14, p = 0.0032). Eight of these SNPs were analyzed for transcriptional regulation using reporter assays where all tested regions exerted regulatory effects on transcription in epithelial and/or immune cell lines, and seven SNPs showed allelic differences. The DNA fragment with the top-associated SNP 27:19,086,778 displayed the highest activity in keratinocytes with 11-fold induction of transcription by the risk allele versus 8-fold by the control allele (pdifference = 0.003), and also mapped close (~3 kb) to an ENCODE skin-specific enhancer region. CONCLUSIONS Our experiments indicate that multiple CAD-associated genetic variants located in cell type-specific enhancers are involved in gene regulation in different cells and tissues. No single causative variant alone, but rather multiple variants combined in a risk haplotype likely contribute to an altered expression of the PKP2 gene, and possibly nearby genes, in immune and epithelial cells, and predispose GSDs to CAD.

Worldwide Diaspora of Aethina tumida (Coleoptera: Nitidulidae), a Nest Parasite of Honey Bees

Native to sub-Saharan Africa, Aethina tumida Murray (Coleoptera: Nitidulidae) is now an invasive pest of honey bee, Apis mellifera L., colonies in Australia and North America. Knowledge about the introduction (s) of this beetle from Africa into and among the current ranges will elucidate pest populations and invasion pathways and contribute to knowledge of how a parasite expands in new populations. We examined genetic variation in adult beetle samples from the United States, Australia, Canada, and Africa by sequencing a 912-base pair region of the mitochondrial DNA cytochrome c oxidase subunit I gene and screening 10 informative microsatellite loci. One Canadian introduction of small hive beetles can be traced to Australia, whereas the second introduction seems to have come from the United States. Beetles now resident in Australia were of a different African origin than were beetles in North America. North American beetles did not show covariance between two mitochondrial haplotypes and their microsatellite frequencies, suggesting that these beetles have a shared source despite having initial genetic structure within their introduced range. Excellent dispersal of beetles, aided in some cases by migratory beekeeping and the bee trade, seems to lead to panmixis in the introduced populations as well as in Africa.

Hybridization between distant lineages increases adaptive variation during a biological invasion: stickleback in Switzerland

The three-spined stickleback is a widespread Holarctic species complex that radiated from the sea into freshwaters after the retreat of the Pleistocene ice sheets. In Switzerland, sticklebacks were absent with the exception of the far northwest, but different introduced populations have expanded to occupy a wide range of habitats since the late 19th century. A well-studied adaptive phenotypic trait in sticklebacks is the number of lateral plates. With few exceptions, freshwater and marine populations in Europe are fixed for either the low plated phenotype or the fully plated phenotype, respectively. Switzerland, in contrast, harbours in close proximity the full range of phenotypic variation known from across the continent. We addressed the phylogeographic origins of Swiss sticklebacks using mitochondrial partial cytochrome b and control region sequences. We found only five different haplotypes but these originated from three distinct European regions, fixed for different plate phenotypes. These lineages occur largely in isolation at opposite ends of Switzerland, but co-occur in a large central part. Across the country, we found a strong correlation between a microsatellite linked to the high plate ectodysplasin allele and the mitochondrial haplotype from a region where the fully plated phenotype is fixed. Phylogenomic and population genomic analysis of 481 polymorphic amplified fragment length polymorphism loci indicate genetic admixture in the central part of the country. The same part of the country also carries elevated within-population phenotypic variation. We conclude that during the recent invasive range expansion of sticklebacks in Switzerland, adaptive and neutral between-population genetic variation was converted into within-population variation, raising the possibility that hybridization between colonizing lineages contributed to the ecological success of sticklebacks in Switzerland.

The long pentraxin PTX3 at the crossroads between innate immunity and tissue remodelling

Innate immunity represents the first line of defence against pathogens and plays key roles in the activation and orientation of the adaptive immune response. The innate immune system comprises both a cellular and a humoral arm. Components of the humoral arm include soluble pattern recognition molecules that recognize pathogen-associated molecular patterns and initiate the immune response in coordination with the cellular arm, therefore acting as functional ancestors of antibodies. Pentraxins are essential constituents of the humoral arm of innate immunity and represent a superfamily of highly conserved acute phase proteins, traditionally classified into short and long pentraxins. Pentraxin 3 (PTX3) is the prototypic member of the long pentraxins subfamily. As opposed to C-reactive protein, whose sequence and regulation have not been conserved during evolution from mouse to man, the evolutionary conservation of sequence, gene organization and regulation of PTX3 has allowed addressing its pathophysiological roles in genetically modified mice, in diverse conditions, ranging from infections to sterile inflammation, angiogenesis and female fertility. Despite this conservation, a number of predominantly non-coding polymorphisms have been identified in the PTX3 gene which, when associated in particular haplotypes, have been shown to be relevant in clinical conditions including infection and fertility. Here we review the studies on PTX3, with emphasis on pathogen recognition, tissue remodelling and crosstalk with other components of the innate immune system.

Estimating the net contribution of interleukin-28B variation to spontaneous hepatitis C virus clearance

The identification of associations between interleukin-28B (IL-28B) variants and the spontaneous clearance of hepatitis C virus (HCV) raises the issues of causality and the net contribution of host genetics to the trait. To estimate more precisely the net effect of IL-28B genetic variation on HCV clearance, we optimized genotyping and compared the host contributions in multiple- and single-source cohorts to control for viral and demographic effects. The analysis included individuals with chronic or spontaneously cleared HCV infections from a multiple-source cohort (n = 389) and a single-source cohort (n = 71). We performed detailed genotyping in the coding region of IL-28B and searched for copy number variations to identify the genetic variant or haplotype carrying the strongest association with viral clearance. This analysis was used to compare the effects of IL-28B variation in the two cohorts. Haplotypes characterized by carriage of the major alleles at IL-28B single-nucleotide polymorphisms (SNPs) were highly overrepresented in individuals with spontaneous clearance versus those with chronic HCV infections (66.1% versus 38.6%, P = 6 × 10(-9) ). The odds ratios for clearance were 2.1 [95% confidence interval (CI) = 1.6-3.0] and 3.9 (95% CI = 1.5-10.2) in the multiple- and single-source cohorts, respectively. Protective haplotypes were in perfect linkage (r(2) = 1.0) with a nonsynonymous coding variant (rs8103142). Copy number variants were not detected. CONCLUSION: We identified IL-28B haplotypes highly predictive of spontaneous HCV clearance. The high linkage disequilibrium between IL-28B SNPs indicates that association studies need to be complemented by functional experiments to identify single causal variants. The point estimate for the genetic effect was higher in the single-source cohort, which was used to effectively control for viral diversity, sex, and coinfections and, therefore, offered a precise estimate of the net host genetic contribution.

Coastal upwelling is linked to temporal genetic variability in the acorn barnacle Balanus glandula

Dispersal and recruitment are central processes that shape the geographic and temporal distributions of populations of marine organisms. However, significant variability in factors such as reproductive output, larval transport, survival, and settlement success can alter the genetic identity of recruits from year to year. We designed a temporal and spatial sampling protocol to test for genetic heterogeneity among adults and recruits from multiple time points along a similar to 400 km stretch of the Oregon (USA) coastline. In total, 2824 adult and recruiting Balanus glandula were sampled between 2001 and 2008 from 9 sites spanning the Oregon coast. Consistent with previous studies, we observed high mitochondrial DNA diversity at the cytochrome oxidase I locus (884 unique haplotypes) and little to no spatial genetic population structure among the 9 sites (Phi(ST) = 0.00026, p = 0.170). However, subtle but significant temporal shifts in genetic composition were observed among year classes (Phi(ST) = 0.00071, p = 0.035), and spatial Phi(ST) varied from year to year. These temporal shifts in genetic structure were correlated with yearly differences in the strength of coastal upwelling (p = 0.002), with greater population structure observed in years with weaker upwelling. Higher levels of barnacle settlement were also observed in years with weaker upwelling (p < 0.001). These data suggest the hypothesis that low upwelling intensity maintains more local larvae close to shore, thereby shaping the genetic structure and settlement rate of recruitment year classes.

Dual origins of dairy cattle farming - evidence from a comprehensive survey of European Y-chromosomal variation

BACKGROUND: Diversity patterns of livestock species are informative to the history of agriculture and indicate uniqueness of breeds as relevant for conservation. So far, most studies on cattle have focused on mitochondrial and autosomal DNA variation. Previous studies of Y-chromosomal variation, with limited breed panels, identified two Bos taurus (taurine) haplogroups (Y1 and Y2; both composed of several haplotypes) and one Bos indicus (indicine/zebu) haplogroup (Y3), as well as a strong phylogeographic structuring of paternal lineages. METHODOLOGY AND PRINCIPAL FINDINGS: Haplogroup data were collected for 2087 animals from 138 breeds. For 111 breeds, these were resolved further by genotyping microsatellites INRA189 (10 alleles) and BM861 (2 alleles). European cattle carry exclusively taurine haplotypes, with the zebu Y-chromosomes having appreciable frequencies in Southwest Asian populations. Y1 is predominant in northern and north-western Europe, but is also observed in several Iberian breeds, as well as in Southwest Asia. A single Y1 haplotype is predominant in north-central Europe and a single Y2 haplotype in central Europe. In contrast, we found both Y1 and Y2 haplotypes in Britain, the Nordic region and Russia, with the highest Y-chromosomal diversity seen in the Iberian Peninsula. CONCLUSIONS: We propose that the homogeneous Y1 and Y2 regions reflect founder effects associated with the development and expansion of two groups of dairy cattle, the pied or red breeds from the North Sea and Baltic coasts and the spotted, yellow or brown breeds from Switzerland, respectively. The present Y1-Y2 contrast in central Europe coincides with historic, linguistic, religious and cultural boundaries.

A shared 336 kb haplotype associated with the belt pattern in three divergent cattle breeds

We recently mapped the belt mutation in Brown Swiss cattle to a 922 kb interval on BTA3. In this study, we analysed two additional cattle breeds with the belted phenotype: Galloway and Dutch Belted (Lakenvelder). By genotyping microsatellites in solid-coloured and belted Galloways, we confirmed that the belt mutation in Galloways is strongly associated with the same chromosomal locus as in Brown Swiss cattle. Subsequently, we analysed 36 SNPs in the belt interval in three breeds. We identified a single belt-associated haplotype for each of the analysed breeds. The three breed-specific belt haplotypes share alleles in four blocks. Three of these blocks comprise only one single or two consecutive markers, while the largest shared haplotype block encompasses nine consecutive SNPs in a 336 kb interval. The large shared haplotype across divergent breeds suggests a common mutation for the belt phenotype in all three breeds. We identified a potential candidate gene within this interval coding for the developmental transcription factor HES6. We re-sequenced the complete HES6 coding sequence in belted and solid-coloured cattle but did not find belt-associated polymorphisms. In conclusion, our data provide strong evidence in favour of a common founder for the belt phenotype in different cattle breeds and have resulted in an improved fine-mapping of the causative mutation.